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We present results on spin dependent potentials from lattice simulations of SU{2) gauge theory. The Coulomb 
like short range part of the central potential is identified as a mixed vector-scalar exchange while the linear long 
range part is pure scalar. 



1. INTRODUCTION 

Potential models have been found to be surpris- 
ingly successful in describing the meson spec- 
trum, even at rather small quark masses, m [Q. 
To derive the semi-relativistic Hamiltonian one 
usually starts from the Bethe-Salpeter amplitude 
of a two-quark bound state with quark separation 
r, takes the static solution (m — > oo) of the ker- 
nel, and expands it in inverse powers of the quark 
mass around this limit. The result can be decom- 
posed into the five relativistic invariants (S,V,T,A 
and P) with r dependent (unknown) coefficient 
functions, the spin dependent potentials (SDPs). 

So far, SDPs have been studied in the mid 
eighties on the lattice. The applicabil- 

ity of the potential approach to heavy quark 
physics can be tested on the lattice by comparing 
the spectrum and wave functions, obtained from 
the semi-relativistic Hamiltonian with SDPs, to 
NRQCD results. In addition we aim at replacing 
phenomenologically adjusted potentials by the 
"correct" QCD potentials. 

The order l/m? Hamiltonian, obtained in the 
manner described above, is a generalisation of the 
familiar Breit-Fermi Hamiltonian. Apart from a 
purely kinetic term, it reads as follows (for equal 
masses), Vsd = Vb + (Vls + ^ss + Kd) /m^, 
with 
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where Li ~ f x pi. Vyd are velocity dependent 
corrections. 

Explicit expectation values have been associ- 
ated to the spin-orbit and spin-spin "potentials" 
V{, V2 and V3, V4, which can be computed on the 
lattice . These potentials are related to scalar 
(S), vector (V) and pseudo-scalar (P) exchange 
contributions in the following way: 
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(r'(r)-P"(r)) (4) 

r 

Viir) = 2S/^Vir) +\/^P{r) . (5) 

In addition, relativistic invariance yields the 
Gromes relation [|, Vo'(r) = V^ir) - V{ir). 

2. LATTICE TECHNIQUES 

From the Cornell form of the central potential, 
Vb(i?) = Vo + KR— e/R, and tree level perturba- 
tion theory one expects 

V; = -K,V^ = ^,V^ = ^,V4 = 8^e<5(i?)(6) 

Since all these potentials (up to V{) are predicted 
to be short ranged, we have chosen small lattice 
spacings to be sensitive to the physically inter- 
esting region: we operate at the values /? — 2.74 
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and /? = 2.96 that correspond to resolutions 
a = .0408(2) fm and a ^ .0183(6) fm, respec- 
tively. The physical scale has been obtained from 
the relation ^Ka = 440 MeV. More than 200 
statistically independent configurations of lattice 
volume 32* have been collected at each coupling. 

In order to reduce statistical fluctuations, tem- 
poral links are integrated analytically wherever 
possible [Q. The SDPs have been defined in a 
way analogous to the computation of local masses 
to exclude corrections proportional to 1/T. 

The lattice correlation functions, C{T), to 
be measured are Wilson loops with two colour 
field insertions (ears) within the temporal trans- 
porters, divided by the corresponding loop with- 
out ears (T denotes the ear-ear temporal dis- 
tance). One ear excites the gluon field between 
the two charges while the second returns the field 
into its ground state. To obtain the SDPs, an in- 
tegration over all possible interaction times (tem- 
poral positions of the ears) has to be performed. 
The minimal distance. A, of an ear to an "end" 
of the Wilson loop, occurring within the integra- 
tion, represents the time the gluon field has to 
decay into its ground state after creation and, 
therefore, governs excited state contaminations. 
The spatial transporters within the Wilson loops 
are smeared Q to suppress such pollutions from 
the beginning, allowing us to reduce A: plateaus 
in the ear-ear correlators have been found from 
A = 2 onwards (in some cases even A = 1). 

In our implementation the positions of the ears 
in respect to the Wilson loop ends are kept fixed 
at A = 2. The interaction time is varied by in- 
creasing the temporal Wilson loop extent. This 
results in smaller statistical fluctuations, com- 
pared to the standard deflnition Q in which the 
temporal extent of the Wilson loop is kept flxed 
and the positions of the ears are varied. 

Previous authors replaced the integration over 
T by a discrete sum. The resulting discretiza- 
tion error adds to other lattice artefacts. We 
take into account that C{T) depends on T in a 
multi-exponential way and approximate the inte- 
gral from T to T -I- 1 by the expression {C{T) — 
C{T + l))/log(C(T)/C(r + 1)), instead. For V( 
and V2 where C{T) is weighted by an additional 
factor T an analogous interpolating formula is 



used. The T — > 00 tail of the correlator is es- 
timated from an exponential fit. An integration 
cut-off Tmax IS choseu such that the statistical er- 
ror exceeds the estimated contribution from the 
asymptotic tail by a factor four. 

The lattice SDPs undergo multiplicative renor- 
malizations in respect to their continuum coun- 
terparts, while Vq does not. A non-perturbative 
renormalization of these quantities in the manner 
suggested in Ref. |^ is performed. In fig. 1 the 
quality of this normalization is checked by com- 
paring the Gromes combination V^' — V( (in units 
of the string tension, K) to the force (sohd curve), 
derived from a fit to the central potential. We find 
good scaling behaviour among the two data sets 
(/? = 2.74 and 2.96) and vahdity of the Gromes 
relation outside of the region of lattice artifacts. 
A more detailed description of the improvement 
tricks and data analysis can be found in Ref. S. 
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Figure 1. Test of the Gromes relation. 
3. RESULTS 

We find the central potential, Vq, to contain a 
very clear linear confining part up to large dis- 
tances which must be of scalar exchange type as 
a vector exchange can grow at most logarithmi- 
cally with r . We confirm the second spin-orbit 
potential V2 to be definitely of short range na- 
ture, leaving little room for a scalar contribution 
to this potential. The first spin-orbit potential, 
V{, is displayed in fig. 2. Our lattice resolution 
enables us to establish an attractive short range 
contribution that can be well fitted to a Coulomb 
(1/i?^) form, in addition to the constant long- 
range term, which is in agreement with the string 
tension, K. 



3 



In principle, the potentials Vb, Vis, and V4 allow 
for a determination of S, V, and P (eqs. 
At this stage, we will assume P to vanish and 
Vi to be pure scalar (which induces the equalities 
V2 — V and Vi = —S). This leads to a prediction 
for V3, according to eq. (||), which can be checked 
in fig. 3: we find reasonable agreement between 
the data sets and the predicted curve, V2/R — 
V2 , the deviations being qualitatively understood 
from tree level lattice perturbation theory. 
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Figure 2. V(, together with a fit curve a/r^ + K. 
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Figure 3. V3 and V^/r - V^' (curve). 

The remaining spin-spin potential, V4, ex- 
hibits oscillatory behaviour as a lattice artifact 
(fig. 4) and can largely be understood as a S- 
contribution, according to V4 = 2V^V2. The tree 
level lattice perturbative expectation 87rc(5L(i?) is 
indicated by crosses in the figure. The normal- 
ization has been obtained from a c/R^ fit to the 
V2 data points. The error ranges without sym- 
bols are obtained by using single lattice gluon ex- 
change with an infra-red protected two-loop run- 



ning coupling 1^. The range corresponds to dif- 
ferent choices of the A-parameter. Apart from the 
dominant (5-like contribution, another very short 
ranged contribution seems to exist. 
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Figure 4. V4 in lattice units. 
4. CONCLUSIONS 



We find a consistent picture when identifying — Vi 
(V2) with the scalar (vector) exchange contribu- 
tion, other exchange types being negligible to or- 
der 1/m?. Apart from the linear large distance 
part, we observe a short distance Coulomb like 
scalar contribution from the first spin-orbit po- 
tential. It appears that the Coulomb part of the 
central potential splits up into a vector/scalar ra- 
tio in between 3/1 and 4/1. The short range po- 
tentials can be qualitatively understood in terms 
of perturbation theory. 
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